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Abstract 

We present a new method for the momentum expansion of Feynman integrals 
with arbitrary masses and any number of loops and external momenta. By using 
the parametric representation we derive a generating function for the coefficients 
of the small momentum expansion of an arbitrary diagram. The method is 
applicable for the expansion w.r.t. all or a subset of external momenta. The 
coefficients of the expansion are obtained by applying a differential operator to a 
given integral with shifted value of the space-time dimension d and the expansion 
momenta set equal to zero. Integrals with changed d are evaluated by using the 
generalized recurrence relations proposed in [Q] . We show how the method works 
for one- and two- loop integrals. It is also illustrated that our method is simpler 
and more efficient than others. 
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1. Introduction 



Accurate theoretical predictions for precision experiments in particle physics are usually 
the result of formidable calculations of radiative corrections. More and more particles 
with different masses have to be taken into account. This requires the evaluation 
of complicated Feynman integrals with many masses and external momenta. Under 
lucky circumstances one can get analytic results for the radiative corrections but more 
frequently one is forced to expand w.r.t. external momenta and/or masses or to find 
other kinds of approximate calculations 0, 0. Frequently even if the result for a 
particular diagram is known analytically [§], the numerical evaluation proceeds via 
approximations which again may be some kind of momentum expansion. It is evident 
that methods of asymptotic expansions of Feynman integrals will play an increasing 
role in the evaluation of physical amplitudes. 

The asymptotic momentum expansions combined with conformal mapping and sub- 
sequent resummation by means of Fade approximations were used to calculate the dia- 
grams in the whole cut plane in the external momentum squared 0-0. As was found 
in these papers the method gives very precise results in a wide range of kinematical 
variables. To extend further this method we need to find a systematic and effective 
method for the momentum expansion of diagrams depending on several external mo- 
menta and masses. For the three-point Green functions this problem was solved in |^ , 
rO[] . As for diagrams with more external legs no reasonable algorithm is known up to 



now. 

The solution of this problem is important not only for the case of small momenta. 
Asymptotic expansions at large momenta or masses also include small momentum 
expansion as an ingredient part ||rT| . 

In the present paper we describe a new approach to the momentum expansion 
which solves the aforementioned problem. The method is closely related to the new, 
generalized recurrence relations for Feynman integrals proposed in 

The paper is organized as follows. In Sect. 2, we present the main ideas of our 
method. First, using the parametric representation, we give a derivation of the gen- 
erating function for the coefficients of the expansion of an arbitrary propagator type 
integral with different masses. Then we derive a similar formula for the diagrams with 
arbitrary number of external legs. 

In Sect. 3 we demonstrate how the method works for one- loop integrals. A general 
formula, in terms of a one-fold integral, for the coefficients of the expansion of scalar 
n-point one-loop diagrams is derived. 

In Sect. 4 several examples of the small momentum expansion of various two- loop 
diagrams are considered. We show how to expand integrals w.r.t. subsets of external 
momenta. The characteristic functions required for expanding propagator type and 
three-point vertex diagrams are given. 



2. General formalism 



The subject of our consideration will be dimensionally regulated scalar Feynman inte- 
grals. An arbitrary scalar L loop integral can be written as 



N 



G^'\{s.}A^l}) = \{Jd%\{Pl^^, (1) 

where 



i=i j=i 



n=l m=l 



Qm are external momenta, {s,} is a set of scalar invariants formed from g^, is the 
number of lines, E is the number of external legs, a; and rj are matrices of incidences 
of the graph with the matrix elements being ±1 or (see, for example, Ref. 

In what follows we will need a parametric representation for G('^)({sj},{m^}). For 
an arbitrary scalar Feynman integral in d dimensional space-time we have 
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where D{a) and Q{{si}, a) are homogeneous polynomials in a of degree L and L -|- 1, 
respectively. They can be represented as sums over trees and two-trees of the graph 
(see, for example, Ref.[ll4|): 



^(«) = E ( n •••«.••), 

over over 
trees chords 

Qm,a) = y: n ( E (4) 

over over over comp. 

2 — trees chords of 2 — tree 

These polynomials are characteristic functions of the topology of the diagram and of 
its subgraphs. Since D and Q will play an important role in the rest of the present 
paper we remind the reader of the definitions of the trees and two-trees for connected 
diagrams. Any connected subdiagram of a diagram G containing all the vertices of G 
but is free of cycles (loops) is called a tree of G. Similarly, a two-tree is defined as 
any subdiagram of G containing all the vertices of the original diagram, but is free of 
cycles, and consisting of exactly two connected components. Finally, a chord of a tree 
(two-tree) is defined as any line not belonging to this tree (two-tree). 

For the practical application of our method the following properties of the functions 
D and Q are important. Namely, the dependence of D and Q on each aiy, can be written 
in the form: 

D{a) = D-u{a)a„ + D^{a), 

Qi{s}, a) = Qu{{s}, a)a^ + Qui{s}, a), (5) 



where D-jj{a), Quds}, a) are the D, Q functions of the diagram obtained by removing 
the hne z/, and D^{a), Qu{{s},<y) are the D, Q functions for the diagram obtained 
by contracting the hne z/ to a point. D^, Du, and Qy do not depend on a^. 
The polynomials -D(a), Q{{si},a) can be easily constructed by means of a computer 



program for any particular integral |T^. Examples of D and Q for some particular 
diagrams will be given in the next sections. 

Let us first consider propagator type integrals. In this case Q is proportional to 
the external momentum squared q^: Q{q^,a) = Q{a)q^. Assuming temporarily that 
all propagators have different masses, from the representation @ we obtain: 



AG('^)(g^K}) = ^-^Q{d)G^'^'\q\{ml}\ (6) 



where Q (d), is the differential operator obtained from Q{a) by replacing aj — > dj = 
d I dm'j . In the above formula it is assumed that after differentiation all masses are set 
to their original value. 

Thus, differentiation w.r.t. momenta of the d dimensional integral is replaced by 
differentiation w.r.t. masses of the d + 2 dimensional one. As we shall show in the next 
sections this replacement in many cases will essentially simplify practical calculations. 
For the l-th derivative we can write 



-1) 



WJ ^""'^^ ^ ^ g'(^)G(^+^')(g^ {ml}). (7) 

Upon setting = we obtain the coefficients of the small momentum expansion of 
the propagator type diagram. 

The general prescription for the small momentum expansion is then the following. 
First, one should assign to all lines different masses and then set the external momenta 
equal to zero. Second, the differential operator Q{d) must be applied as much as needed 
and then the masses must be set to the original one's. Third, the integration by parts 
recurrence relations in c? + 2/ dimensions should be used, reducing all the integrals to 
the set of master integrals. Fourth, the set of all master integrals in + 2Z dimensions 
must be expressed in terms of d dimensional master integrals using the generalized 
recurrence relations described in [|l| . In some cases it is advantageous to use recurrence 
relations which simultaneously reduce the value of d and the sum of exponents of the 
propagators in the integral. Rather frequently master integrals are just combinations 
of Euler's F functions, and transforming the latter from d + 2l to d dimensions is quite 
simple. 

In practical calculations no real differentiation is needed. It is enough to compute 
Q'(a;), multiply it by G^'^+^')(0, {m^}) and make the substitutions 



[A;2 - mf + ie]" F(n) [k^ - mf + ie]"+J' ' 



An important feature of our method is that one does not need to care about ex- 
ternal momenta, avoiding many complications intrinsic to other methods based on the 
differentiation w.r.t. external momenta. From the very beginning, one is dealing with 



simpler objects; only bubble diagrams need to be considered. In all cases where the 
bubble diagrams are known, one immediately gets an explicit generating function for 
the coefficients. 

Unfortunately the polynomial Q usually has many terms, and it is not so easy 
to exponentiate it. For the three- loop nonplanar master integral occurring in the 
QED photon propagator, Q is a fourth order polynomial with 45 terms. In this case 
using FORM |jl6|] running on a DEC 3000 it took about 46 hours to generate the 12-th 



coefficient in the small momentum expansion in terms of bubble integrals. Nevertheless, 
the number of generated terms, and as a consequence the execution time, is smaller 
than what one gets using the direct small momentum expansion methods. Using the 
method proposed recently in fl^ for solving recurrence relations for three-loop vacuum 



integrals one may expect an improvement of the efficiency in the complete calculations 
of the coefficients. Our method is definitely the most optimal for the case of all different 
masses. In case when some masses in the integral are equal the calculation can be 
simplified by taking into account this information from the very beginning. 

The algorithm described above can be generalized to the small momentum expan- 
sion of integrals with an arbitrary number of external momenta. For integrals depend- 
ing on several external momenta the polynomial Q is a sum over a set of n = E{E — l)/2 
independent scalar invariants sj: 

n 

Qia) = J2sMa), (8) 

where E- is the number of external legs of the diagram corresponding to a given integral. 
In this case the analogue of (|^), valid for propagators, reads: 

^G('\{s,}, K}) = ^^^^^Q5(a)G('^+2')({.,}, K}). (9) 

To find the coefficient of the M-th order in s''i . . . sj^ one should apply the product 
of operators Q''i (d) . . . {d) to the integral taken in d + 2M dimensions. If we perform 
an expansion w.r.t. all scalar invariants again the coefficients of the expansion will be 
generated from the bubble integrals. In this case the general formula for the small 
momentum expansion of an arbitrary scalar integral will be: 



M=0 " 
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X QHd)...Ql-{d)G('-^'^^\{0},{m^})-^j...-^. (10) 

{h + ... + ln=M} '1- 

The expansion can also be performed w.r.t. some subset of scalar invariants which 
are proportional to the momenta considered to be small. In this case the operators 
Qj [d) must be applied to the given integral where only that subset of the momenta is 
taken to be zero. The corresponding formula is almost identical to (0) and may be 
readily written as: 



TT-ML 

M=0 " 



X Y: Qiid) ■ ■ ■ Ql-{d)G^'^'''\{0}, {S^}, K})^ ...fj, (11) 

{h + ...+l„=M} '1- 

where {Sp} is a set of scalar invariants assumed to be fixed. 

Equations (|l^) and (11) can be used not only for the small momentum expansion. 
The momentum expansion considered is also an essential ingredient of other kind of 
asymptotic expansions 0, |Tl|, [|1^]. 



3. Expansion of one- loop integrals 

Now we shall show how the method works by considering several simple examples. We 
start with the one-loop integral: 



pa p/3 
d/2j ki,mi ki-q,m2 



OO T 



/=0 ''• 
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In this case Q{q^-, o) = a\a2 q^- The /-th coefficient of the small momentum expansion 
will be 

Ql Ql ^ d'^+^^kl 



Upon differentiation we obtain: 



7712 ■ 



T{l + a)Til + P) I 



T{a) T{f3) 



HIT 



d/2+«j fci,"!!-* fci,m2- 



The bubble integral can be expressed in terms of hypergeometric function 

2F1 



d^p. p, _ (-ir^r(a + /3- 

(i/2j fci,mi-' fci,m2 



(m?)"-i(m2)/3 T{a + p) 
For simplicity we will take mf = m| = m^. In this case 



^^^0 r{l + a)r{l + (3) T{a + (3 + 1-1) 



a + (3; ml 



r(a) r(/3) r(a + /5 + 2/)(m2)'+°+/3-'^/2- 



and therefore 



-l)"+^r(a + /5- 



[m 



2)"+/5-lr(« + /?) 



a,P,a + P - _(f 



a+P a+/3+l 
2 ' 2 ' 
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(13) 



(14) 



(15) 



(16) 



(17) 



in agreement with Ref.[19 



To compare our method with it is instructive to calculate the coefficient Ji at 
= = m^, a = P = 1 using both techniques. In our approach Ji is just one term: 



1 



dml dml J [zvr^/^+i] [kj — mf + ie){kj — + ie) 

d'^+^ki 1 

(A;f -m2 + ie)4' 



'tj=m2=m2 
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The same coefficient in the method pO[ is given by 



1 



2d^q^^q^^ J [iTr"'/^] ((/ci — g)^ — + ze)(/cf — + ze) 
1 /• d-^A;! [Akf-d{kl-m^)] 



q=0 



(19) 



One can observe at ffist, that differentiating w.r.t. masses is easier than w.r.t. external 
momenta. Secondly, after differentiation the number of terms for Ji in ( pJl) is lager than 
in our approach. The number of terms for the higher coefficients Ii in the approach 
grows exponentially. A very similar situation we find for integrals with more loops. 
For d dimensional one-loop integrals corresponding to diagrams with E legs the 
number of terms in Q is equal to the number of independent scalar invariants n = 
E{E — l)/2. In this case Q may be represented as: 



Q{{s},a) = Y,Qj{(^)sj, 



(20) 



with Qj{<y) being monomials in a. In the small momentum expansion of the one- loop 
integrals: 



p^E-l pVE 
ki-qi,mi ■ ■ ■ ■'■ ki-qE-i,mE-i ki.rriE 



M=0 {Zi+.„+z„=Af} 

the coefficients of the expansion will be 



hi 



(21) 



a,..i„ = :^QHd) . . . Qi:id) I d'^^'^^k, p^,i^^ . . . p-^^. 



(22) 



Performing the angular integration we obtain a compact one-fold integral representa- 
tion for the coefficients of the momentum expansion of an arbitrary one-loop integral: 



r(f + M) 



(x+mf)'"! . . . (x+m|;) 



(23) 



The previous integral may readily be evaluated in terms of the Lauricella function 
F^~^^ (see formula 7.1.1.5 in Ref. |^) in agreement with Ref. In the case when 



all masses are equal or some masses are zero while the others are equal the integral can 
be evaluated in terms of Euler's F function. 

For illustration, we consider the expansion of the three-point vertex function w.r.t. 
external momenta Q'l, ^3 = (?i + ^2)^: 
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ki—qi,mi ki,m2 ki+q2,m3 / ^ hhl; 



k\ y y ' 



(24) 



The function Q corresponding to IuiU2U3 reads: 

Q{{s}, a) = aia2ql + a2azql + ai^sgg. (25) 
The coefficients of the expansion in q^, g| may readily be found from (|2^ and (|23|): 
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where (z/); = F(z/ + Z)/F(z/) is the Pochhammer symbol. For the case of equal masses 
ml = = ml = m^, the coefficients are just products of Euler's F functions: 



^y,+U2+U3 i^^)h+l3il'2)h+l2il^3)l2+l3^ ('^l + ^^2 + ^^3 + ^1 + ^2 + - j) ^^^^ 



[m 



2y,+u2+U3+h+l2+l3—2 Y{u, + Z/2 + Z/3 + 2(/i + k + h)) 



Let us consider now an integral with one small momentum, and the other ones 
arbitrary. This kind of kinematics occurs, for instance, in the evaluation of the moments 



of the deep inelastic structure functions p3[. The one-loop scalar diagram which is of 
interest in this case is shown in Fig.l. 



q 


1 


q 


4 




2 


p 


3 


p 



Fig.l Box diagram depending on two momenta 



The contribution to the moments of the structure functions from the integral related 
to the diagram in Fig.l will be determined by the coefficients of the expansion w.r.t. 
the scalar product (pq) with the assumption p"^ = 0: 

^(^'' ^) = y T-^^PLrnPk,+q,rnPk,-p,rn = EiP^TTn- (28) 



n=0 



From the function Q of the diagram we derive the differential operator that generates 



the coefficients T„: 



Qp,(9) = 2did^. (29) 
Introducing auxihary masses, settings = 0, applying Qi,Jd) to the diagram and setting 



m 



after differentiation, we get: 



T 



Jd+2n I. 

cyn \ I " pn+3 pn+1 

' r.d/2+n] ki,m-^ki+q,m 



2"n! r (n + 4 - f 



2]n+4-d/2 



■ 3-^2 



n + 3,n + l,n + 4 — 



2 ' 



(30) 



The integral in (|30D was evaluated with the aid of ([17|) . A different method |2^ to 
calculate moments of this kind of diagram, consists in applying a differential projection 
operator built up from the momenta p and q. From the above comparison made for 
the integral ll^p{q^ ^m^^rn^^ the analysis given in the next section for the two-loop 
diagrams and the remark in Sect. 2 concerning three-loop calculations, we may expect 
our method to be more efficient than the method of projections. 



4. Expansion of two-loop integrals 

Small momentum expansions of various kinds of two-loop diagrams have been consid- 
ered in ||20|. For the generic two-loop propagator type diagram 

T |'^2^ ^ f [ d'^ki d'^k2 ^3 , > 

Q{q'^,a) and D{a) are, respectively, the third and the second order polynomials in a: 

Q{q^, a) = Q{a)q^= [(ai + a2)(a;3 + a4)a5 + aia2(a3 + a4) + a3a4(ai + a2)]Q'^,(32) 
D{a) = a5(ai + ^2 + as + 0:4) + («! + a3)(a2 + a^)- (33) 

If, at q^ = 0, JviU2U3V4,u^X'i^) does not have singularities then it can be expanded into 
the Taylor series: 

00 

^...2.3.4.5(?')=E^/ (34) 

The differential operator generating J/ will be a four-fold sum: 



Q\d)=T.T.m: n ■ : ; ; ^ \ : , , n ^, • (35) 



jj^^^^^ - j)' (i-O' {3 -ry. i\ r! s\ {I ~ i - s) 
Applying this operator to the diagram taken at g = we find: 



P P ^d+2lh. Jd+2lh. 

y- 1 p^i+i+s pV2+l+s pu-i+l-j pvi+r+l-j py^+i-i /or.\ 

I I ]^^d/2+l'Y ki,mi ^k2,m2 ^ ki,mz ^k2,mA ^ ki-k2,mz' W^/ 



If z/j = 1, . . . , 4 ) are integer, performing partial fraction decomposition like 

1 



{kf -m\ + ieY^ {kj - m§ + ie^^ 



j=0 

us,' 



j\ (z/3 - 1)! {ml - mlY'''+^{kl - mj + ie)"^-^ 



we can express J/ as a multiple sum each term of which is proportional to a bubble 
integral with only three denominators raised to some powers. These integrals can 
be expressed in terms of four Appell's function F4 which in turn are double 
series. Representation (|36D may be helpful for numerical calculations using the multiple 
precision program ||25|| . 



j{d+2M) ^ " ^1" ^2 pv-i (nc 

v\V2V3 j j L-„d/2+M12 k-i,mi ki-k2,m2 k2,m3^ W'- 



It took 5 min on a PC Pentium 90, using FORM |]I6[, to generate in terms of 
bubble integrals the 30-th coefficient of the small momentum expansion for the two- 
loop master diagram with all masses different. Again, as we observed in the one- loop 
case, the number of generated terms in our approach is smaller than in the direct small 
momentum expansion. For instance, the 1-st, 2-nd, 3-rd, 4-th and 5-th coefficients of the 
Taylor expansion of the two-loop diagram with all masses equal have 3,6,10,15 and 21 
terms, respectively, in our approach and correspondingly 5,18,45,97 and 182 terms in 
the other methods [|], To reduce the bubble integrals 

encountered in the evaluation of the M-th coefficient of the expansion, to the ci = 4 — 2£: 
dimensional ones, it was enough to use only one recurrence relation which is given in 

^^2^/3^4^.?+! .3+1 + '^^i^JuXi U2U, -{d- 2z^i)41.3 = 0- (39) 
Needless to say that in our method the coefficients were expressed in terms of integrals 
with the sum of exponents of scalar propagators larger than in the other methods, 
i.e. more complicated integrals. Despite of this fact, the sum of those integrals were 
calculated faster. For example, the 6-th coefficient of the aforementioned diagram was 
calculated two times faster. For higher coefficients, our method becomes even more 
effective comparing with methods known up to now. Significant improvement in the 
efficiency of the method may be achieved by implementing it as a numerical procedure 
proposed in using the multiple precision program . 

In some simple cases the method proposed may be helpful in finding analytic results. 
For example, the two-loop diagram with three denominators: 

T^''\q\mlmlml) = jj d%d% ^,^^^^,,^,^^,,^3' (40) 

has a rather simple Q function: Q{q'^,a) = a\a20fi(^ ■ The coefficients of the small 
momentum expansion of the diagram: 

00 rp 

T^'\q\mlmlml)=Y.Mq^Y (41) 

1=0 



can be written as: 



(42) 



Using formula (4.3) from Ref. we obtain: 

" r(z.i)r(z.2)r(/.3)r(f + /) 

X |r(^ - Z/i)r(^ - U2)T{UI + 1^2 + 1- ^)T{U, + Z/2 + 1/3 + / - rf) 

r^/ 77 7^ , d ^ d . . 
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+ yi-^'T{^ - z/i)r(z/2 - ^)T{v, + /)r(i.i + ^^3 + / - ^) 

(i d d 

X ^4(1/1 + /, z/i + z/3 + / - -; z/i + 1 - -, - + 1 - i^2|a;, y) 

+ a;f-T(z/i - ^)r(^ - z/2)r(z/2 + /)r(z/2 + ^3 + / - ^) 

7^ / 7 7 d d ^ ^ d . . 

X ^4(1^2 + ^^2 + ^^3 + f - 2 i 2 ^ '^i' + 1 - 2 I ^' ^) 
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r(z/i--)r(/--)r(z/3 + /)r(- + 



d d d 
X F4(z/3 + /, - + /; - + 1 - z/i, - + 1 - z/2 I a;, y) 



2 / 2 
1712 /m^ 



(43) 
By 



where F4 is Appell's hypergeometric function, x = mf/m^ and y — /;t2/"t3- 
inserting ( ^^ into (^) and using the series representation for F4 we reproduced the 
exphcit resuh given in [Q. 

At the two-loop level, there are two topologically different three-point vertex dia- 
grams shown in Figs. 2a and 2b, where each line corresponds to a scalar propagator 
with an arbitrary exponent. 
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Fig. 2 Two — loop vertex diagrams 



The Q{{s}, a) and D{a) polynomials for the planar diagram Fig.2a are: 
Qpi{{s}, a) = {aias + aia^ + 0203 + a3a6)«59i 

+ [aia2a5 + aQ{a2 + 04) ("i + "3) + cna3(tt2 + "4) + 0204(01 + 03)] g|,(44) 

Dpi{a) = {ai + a2)(a3 + ^4 + as) + ce^^ai + 0^2 + 03 + "4 + as)- (45) 
For the nonplanar diagram Fig.2c Q{a) and D{a) are: 

Qnpii{s}, a) = [aittsOe + «2«3 «4 + «3«s(«i + a2 + 04 + ae)] 
+ [aiasa^ + q;2Q;s«6 + 0:406(01 + 02 + 03 + 05)] ^2 
+ [02O3O6 + 01O4O5 + 0102(03 + 04 + 05 + oe)] g|, (46) 

Dnpi{a) = (oi + 02)(o3 + 04 + 05 + oe) + (03 + 05)(o4 + oe). (47) 

Any two-loop scalar diagram at zero momenta G^^\{0} , {ml}) can be written, 
using partial fraction decomposition (^), as a multifold sum over bubble integrals 
with three denominators only. As we already mentioned those integrals for arbitrary 
space-time dimension and arbitrary exponents of propagators are combinations of four 
Appell functions F4. Therefore, using (|TUp, the coefficients of the Taylor expansion of 
any two-loop integral can be written in a closed form as a multifold sum over functions 
F4 times product of Euler's F functions. Multiple sums will be produced by (d) and 
partial fraction decomposition of the propagators with the same momenta and different 
masses. 

The representation of the coefficients in terms of F4 functions may not be of great 
use in general. For the cases when some masses are equal or vanish one can get 
simpler representations. If, for example, two masses in the vacuum integral with three 
denominators are equal, then the F4 functions are reduced to the somewhat simpler 
4F3 functions |^D[. If one mass is zero then F4 reduce to 2-^1 functions. 

As an example, we consider the small momentum expansion w.r.t. to the subset of 
scalar invariants g|, 52^3 for the diagram Fig. 2c: 

ff d'^kid'^k2 

. ff^lMllatoS. (48) 

By using (^, the polynomial Q{a) for this diagram 

Q{{(li(lj}i ") = 2(^2^3)010203 + 0102(03 + 04)^2 + g303[oi02 + 04(01 + 02)], (49) 

can be obtained from (^ID or (H) by setting 03 = 04 = and then replacing —^03, 
Oe — > 04. The coefficients Cir{q2) can be represented as: 

C.M^^i-ry-r.Y ^^_%,_^f (50) 



^d+2r+2l J^^^d+2r+2l ^ 



2 p p p p 

ki—q2,Tni-^ fci— fc2,m2-' fe2,m3 k2,m4- 



By differentiating and then performing partial fraction decomposition, the coefficients 
Ciriqi) will expressed as combinations of integrals: 

„ „ Jd+2r+2l u Jd+2r+2l t. 
j{d+2r+2l) ^ ^ ^1 " lip!^! P'^a pi^a 

'^uiU2U3 J J |'^^(i/2+Z+rj2 -'^fci-g2,nii fci-fc2,m2 fc2,nij5 W-"-/ 

where rrij corresponds to m^, or m^. To reduce d + 2r + 21 dimensional integrals to 
d dimensional the generalized recurrence relations proposed in [|I| must be used. Just 
for illustration we derive several relations of this kind. In particular one may get the 
relation: 

Vl(l2t.jf\^+l U2U3 + {d-Ui-U2~ Vz)JI%u^ 

- mluiji%^ ^^^^ - m^z/sJ^fi^+i - ^l^34fu2^,+i = 0, (52) 

where 

Following the procedure described in [|l| the integral J^^j^^^^^^ can be expressed in terms 
of a scalar integral with the space-time dimension + 2: 

By inserting into (^21) we obtain: 



+ mlviffXi V2VZ + ^\^ijf^v2^\ + mlu34% ^3+1. (55) 

From the equation connecting d and d — 2 dimensional integrals given in which in 
our case reads 

= (^1^2 + + 9293)41.3 (56) 

another relation follows: 

j(d-2.) j{d) j{d) j{d) _ ^ /KvN 

'^1^11^21^3 ~ ^1^2Jui + l U2 + 1 U-j, ~ ^l^S'Jui+l U2Ua + l — ^l^'iJuiU2 + l vz + 1 — U. {0 1} 

By using Eqs. (|55| ) and (|57D, several others obtained in a similar manner and relations 
derived from the method of integration by parts |2^ , one can reduce the integrals ( pl|) 
to the set of c? = 4 — 2e dimensional master integrals. 

As a final example we consider the small momentum expansion of the two-loop 
scalar vertex diagram occurring in the process H — » 27: 

d'^ki d'^k2 

]^^d/2']2 ^ki+qi,mPki+q2,TnPk2+qi,mPk2+q2,mPki-k2,oPk2,rn') (^8) 

with Qi = qi — ^- ^his kinematics the method of integration by parts gives: 

^ki+qi,mJ^^k\+q2,m^k2+qi,mJ^^k2,m [-< A;2+g2,»n ~ fci-fe2,oJ • I'^Jj 



(^d — A) J J [m'^l'^Y '^2+<?i,m-rfc2,m 



The first term is just the product of factorized one-loop integrals which are known in 
terms of hypergeometric functions |jl9|. By using our method we expand the second 
integral in ( ^91 ) w.r.t. g| = (gi — ^2)^ • The Q function can be obtain from (44) by 
setting q'l = q'l = and 04 = 0: 



Q{qj, a) = Q?Xa)ql = «2[ai(a3 + "5) + a6(«i + "s)]?!- 
The differential operator for the l-th term of the expansion will be: 

^1 I i-j j d''~''~^^d'~'^~^d'^d-' 
Qiid) = ^ 4 ggg 0--A:V.! kl- 

By applying this operator to the integral 



(60) 



(61) 



' P P 



kl ,mi ^k2,m2 ^ki ,1x13 ^ki,m5 ^k\—k2,me 1 



P P 



setting mi = m2 = = m, ttlq = and adding the contribution from the factorized 
integrals we obtain the series expansion of / in g| = (gi — ^2)^: 



m 



,12-2d 



Ell 



1=0 



m? 



(62) 



where 



2(/ + l)! 
(rf-4) 



^ r(/-j + 3-f)r(j + 3-f; 

(/ + 1)! 4'+2 A. r(/ + |)r(j + 1) (j + 1) 



T{l + Q~d) 

r(f + 



(63) 



' j! r(f + /-j-i) ^ (/-fc)! r(fc + 3-f)r(/ + j-fc + 4-f) 

""^0 r(/ + j + 7-rf) (z-^-)! (^--A;)! (fc + l)! (2/-fc + 2)! 



For (i = 4 this formula confirms the coefficients given in 



5. Conclusions 

We presented a new systematic method for the momentum expansion of arbitrary scalar 
multiloop Feynman diagrams. The method can be used for the small momentum 
expansion and as an ingredient for the large momentum expansion or other kind of 
asymptotic expansions. An important feature of our method is that it can be used in 
cases when only some momenta are small and the others arbitrary. This extends the 
applicability of our method to a wider class of physical problems. Of special interest 
will be the application of our method to the evaluation of the moments of structure 
functions in deep inelastic scattering. 
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